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\ We study the energy and momentum of an isolated system in the tetrad theory 
of gravitation, starting from the most general Lagrangian quadratic in torsion, which 
>• ' involves four unknown parameters. When applied to the static spherically symmetric 
case, the parallel vector fields take a diagonal form, and the field equation has an exact 
O solution. We analyze the linearized field equation in vacuum at distances far from the 
^ . isolated system without assuming any symmetry property of the system. The linearized 
■ equation is a set of coupled equations for a symmetric and skew-symmetric tensor fields, 
^ but it is possible to solve it up to 0(l/r) for the stationary case. It is found that the 
O . general solution contains two constants, one being the gravitational mass of the source 
I . and the other a constant vector B^. The total energy is calculated from this solution 
tijr)' and is found to be equal to the gravitational mass of the source. We also calculate the 
j>I _ spatial momentum and find that its value coincides with the constant vector 3^- The 
linearized field equation in vacuum, which is valid at distances far from the source, 
does not give any information about whether the constant vector Ba is vanishing or 
not. For a weakly gravitating source for which the field is weak everywhere, we find 
that the constant vector Ba vanishes. 



§1. Introduction 



Einstein concerned himself with the problem of energy and momentum for a sys- 
tem of matter plus gravitational field shortly after he proposed general relativity.^) He 
introduced the well-known expression for the energy-momentum complex 0^^ , which is 
referred to as a canonical expression and satisfies ordinary conservation law as a con- 
sequence of the gravitational field equation. For a closed system and for a restricted 
class of coordinates the quantities obtained from the complex 6^^ by integrating 
over spatial coordinates satisfy the following propertiesi^^'^^ (1) They are constant in 
time, (2) transform as a 4-vector under linear coordinate transformations, and (3) are 
invariant under arbitrary transformations of spatial coordinates which tend to the iden- 
tity transformation at infinity. Furthermore (4) the gravitational mass of the system 
is equal to the total energy. It has recently been shown that a characteristic property 
of energy thus defined is its positivity, reflecting stability of the physical system. 

M0ller revived^) the issue of energy and momentum in general relativity, and 
required^) that any energy-momentum complex r^*^ must satisfy the following prop- 
erties: (A) It must be an affinc tensor density which satisfies conservation law, (B) 
for an isolated system the quantities are constant in time and transform as the 
covariant components of a 4-vector under linear coordinate transformations, and (C) 
the superpotential U^^^ = —U^'^'^ transforms as a tensor density of rank 3 under the 
group of spacetime transformations. The second property requires that t^'' around an 
isolated system should decrease faster than for large values of the radial coordi- 
nate r. The last property is required in order that the four-momentum for any isolated 
system be a 4-vector under arbitrary spacetime transformations. 

The canonical energy-momentum complex dfj^ does not satisfy (C). An energy- 
momentum complex that satisfies properties (A) and (C) was constructed,^'*'^'' but it 
was found later^'^) that the property (B) is not satisfied. Lessner^^^ discussed that 
although this result is inadequate from the viewpoint of special relativity, it may be 
adequate from the viewpoint of general relativity. 

It is not possible to satisfy all the above requirements if the gravitational field is 
described by the metric tensor alone. In a series of papers, ^°-''^^)'^^-' therefore, M0ller 
was led to the tetrad description of gravitation, and constructed a formal form of 
energy-momentum complex that satisfies all the requirements. The metric tensor is 
uniquely fixed by the tetrad field, but the reverse is not true, since the tetrad has six 
extra degrees of freedom. In the tetrad formulation of general relativity, the tetrad field 
is allowed to undergo local Lorcntz transformations with six arbitrary functions. The 
energy-momentum complex is not a tensor and changes its form under such transfor- 
mations. Therefore, unless one can find a good physical argument for fixing the tetrad 
throughout the system, one cannot speak about the energy distribution inside the 
system. The total energy-momentum obtained by the complex, however, is invariant 
under local Lorentz transformations with appropriate boundary conditions. 

M0ller also suggested another possibility that the tetrad field is uniquely fixed 

by means of six supplementary conditions. In this case the underlying spacetime 
possesses absolute parallelism,^"^) with the tetrad field playing the role of the parallel 
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vector fields, which are allowed to undergo only global Lorentz transformations. The 
Lagrangian formulation of this tetrad theory of gravitation was first given by Pellegrini 
and Plebanski.^^) Hayashi and Nakano^^^ independently formulated the same gravita- 
tional theory as a gauge theory of spacetime translation group. The Lagrangian was 
assumed to be given by a sum of quadratic invariants of the torsion tensor, which 
is expressed by first-order derivatives of the parallel vector fields. Its most general 
expression then involves four unknown parameters to be determined by experiment, 
which we denote here as ai, 02, 03 and a^. The last parameter 04 is associated with a 
parity-violating term. At first it was required that for the weak field case the gravita- 
tional field equation should reproduce the linearized Einstein equation. This restricted 
the parameters as ai -|- 02 = and 04 = 0.^^-''^''^ M0ller^'^) also suggested a possible 
generalization of the gravitational Lagrangian by including homogeneous functions of 
the torsion tensor of degree 4 or higher. 

Hayashi and Shirafuji^^-* studied the geometrical and observational basis of the 
tetrad theory of gravitation^^ assuming the Lagrangian to be invariant under parity 
operation, involving three unknown parameters ai, 02 and 03. Two of these parameters 
Oi and a2, were determined by comparison with solar-system experiments, while an 
upper bound was estimated for the 03 .Q^ It was found that the numerical value of 
ai + 02 should be very small, consistent with being zero. 

In the tetrad theory of gravitation, as far as we know, the total energy of an 
isolated system has been calculated only for spherically symmetric case. In the case 
di + 0,2 = = Mikhail et al.^^-* found a static, spherically symmetric solution (6^^) 
in Cartesian coordinates with (b^^^a) ~ l/v^ ~ (^"0) for r— > 00, and showed that the 
total energy does not coincide with the gravitational mass. This result was extended to 
a wider class of solutions with spherical symmetry. An explicit expression was given 
for all the stationary, asymptotically flat solutions with spherical symmetry, which 
were then classified according to the asymptotic behavior of the components of {b°'o) 
and (b^^^a)- It was found that the equality of the gravitational and inertial masses 
holds only when (&"o) and (b^'^^a) tend to zero faster than l/i/r. 

Calculation of the energy of an isolated system was extended to the generic case 
(ai + a2)(ai —Aas/9) 7^ under the assumption of spherical symmetry."^^^ It was shown 
that linear approximation can be applied to the field equation at large spatial distances 
from the source, and that the calculated energy is equal to the gravitational mass of 
the source. When (ai — 403/9) = 0, however, the problem of energy of an isolated 
system is not yet understood well, even for the spherically symmetric case, since it is 
not always possible to determine the asymptotic behavior of the components by 
means of linearized field equation. 

It is the purpose of this paper to calculate the total energy and spatial momentum 
of an isolated system without assuming spherical symmetry. For this purpose we solve 
the linearized field equation at far distances up to order 0(l/r). The general solution 

*They coined the name "new general relativity", because Einstein^^' was the first to introduce the 
notion of absolute parallelism into physics after he had constructed general relativity. 

For macroscopic matter, Nitsch and Hehl^"^ proposed a tetrad theory of gravitation as the 
translational guage limit of Poincare gauge theory. Their choice of the parameters corresponds to 
ai = —1/3, 02 = 1/3, 03 = 3/2 and = in our notation. 



3 



has two constants. One of this is related to the gravitational mass of the source, and 
the other is a constant vector. We then show that the energy is always equal to the 
gravitational mass of the isolated system. In this sense the equivalence principle is 
satisfied in the tetrad theory of gravitation. It should be noted that our conclusion is 
based on the assumption that weak field approximation can be applied at far distances 
from the source. In the spherically symmetric case mentioned above, however, solutions 
are known for which the linear approximation does not work even far from the source 
because the components, (fo'^o) and (b^^^a), behave like I/a/t. 

In §2 we discuss the most general Lagrangian with a parity- violating term and apply 
its field equation to the static, diagonal parallel vector fields with spherical symmetry. 
We find that the exact solution for the parity-conserving case^^-* satisfies the field 
equation. In §3 we construct the linearized form of the field equation in vacuum 
at distances far from an isolated system. Assuming that the system is stationary, 
we obtain the general solution of the linearized field equation in vacuum. In §4 we 
derive the superpotential from the general Lagrangian and calculate its components 
necessary for computing the total energy and momentum of the system, using the 
solution obtained in 53. The final section is devoted to conclusion and discussion. 



§2. Basic Lagrangian 

In a spacetime with absolute parallelism the parallel vector fields (b^^) [f define the 
nonsymmetric connection 

= bk'b\^, (1) 

with b^^^u = dub'^f^, from which the torsion tensor is given by 

T /^i/ = r r ufi ~ bk (b b u,^- (2) 

The curvature tensor defined by T^p,i, is identically vanishing, however. Here Latin 
indices are raised or lowered by the Minkowski metric riij=ri^^ =diag(— 1, +1, +1). 
The metric tensor is given by the parallel vector fields as 

g^^u = bk^bK. (3) 



Assuming invariance under a) the group of general coordinate transformations, and 
b) the group of global Lorentz transformations, we write the most general gravitational 
Lagrangian density quadratic in the torsion tensor asQ 



G 



(4) 



*Latin indices {i, j, k, ■ ■ ■) designate the vector number, which runs from (0) to (3), while Greek 
indices (/i, v, p,- ■ ■) designate the world-vector components running from to 3. The spatial part of 
Latin indices is denoted by (a, 6, c, • • •), while that of Greek indices by (a, /3, 7, • • •). 

^ Throughout this paper we use the relativistic units, c = G = 1. The Einstein constant k is then 



equal to Stt. We will denote the symmetric part by ( 
the antisymmetric part by the square bracket [ ], Aj^^j 



for example, ^(;^^) 



{1/2){A^,^ + A^^) and 



4 



where ai, a2, as and are dimensionless parameters of the theory, and t^^A, Vfi and 
are the three irreducible components of the torsion tensor. ^^^'0 

By applying the variational principle to the Lagrangian (4), we obtain the field 
equation: 

with 

I"^ = 2k[DxF''^^ + vxF"^^ + H"'' - ^g^^Lc], (6) 

where 

1 n r 

2 dbK,, ^ ' ^'^ 
= TP^Tp^" - ^T^^^F^p, = H"^, (8) 

J.,. = (9) 

Here Lq = £^G/y/~9, and Cm denotes the Lagrangian density of material fields, of 
which the energy-momentum tensor T^^" is nonsymmetric in general. 

In static, spherically symmetric spacetime, the parallel vector fields take a diagonal 
form, and the field equation (5) can be exactly solved. The exact solution so obtained 
is the same as that obtained by assuming the invariance under parity operation. ^^-^ This 
implies that the parameter 04 has no effect when we use diagonal parallel vector fields 
having spherical symmetry. 

In order to reproduce the correct Newtonian limit, the parameters ai and 02 should 
satisfy the condition 

ai + 4a2 + 9aia2 = 0, (10) 
called the Newtonian approximation condition, ^^-^ which can be solved to give 

"^ = -3(1^' "^ = 3(1^' 

in terms of a dimensionless parameter e. Comparison with solar-system experiments 
indicates that |e| must be very small. 



§3. Solution at far distances 



In order to calculate the energy and momentum of an isolated system confined in 
a finite spatial region, we focus our attention to the solution at far distances from 
the source without assuming spherical symmetry. In such spatial regions far from the 
source, the gravitational field is weak and matter fields are not present: Thus, we are 

The dimensionless parameters KUi of Ref. 18) are here denoted by for convenience. 
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allowed to treat the linearized field equation in vacuum. For simplicity we assume that 
the whole spacetime is covered with a single coordinate system {x^; +00 > > —00} 
with the origin being located somewhere inside the finite system. In spatial regions far 
from the source, the parallel vector fields can be represented as 

h\{x) = 5\ + a\, |aM«l, (12) 
hk^{x) = 5k^ + Ckr |cfe1«l, (13) 

and we assume that all the quadratic and higher-order terms of au^ (or c^^) can be 
neglected in the field equation. Accordingly, we do not distinguish Latin indices from 
Greek indices for ^ and Ck^: We use Greek indices which are now raised and lowered 
by the Minkowski metric ■q'^^ and 77^,^. ^Prom the relation ^jbu" — ^tx^ , it follows that 

+ c^i. = (14) 

with a^^ = r)u\ci^n and c^i, — rj^xCn^- We decompose a^,^ into symmetric and antisym- 
metric parts, 

with h^i, — hi,fj, and A^^, = —Ay^. The components of the metric tensor are written as 

Qiiv = Tlt.y + h^y. (16) 

The antisymmetric part has no contribution to the spacetime metric, implying that it 
is associated with the intrinsic spin-1/2 fundamental particles. 

Keeping only the linear terms and putting T^j^ = on the right-hand side, we see 
that the symmetric part and the skew-symmetric part of the field equation (5) take 
the form 



— n V - 2^"i + «2)(^^:/n/i - h^fj,^) + (y - a2)Vnuh ,pa- 

_ 2 

(2ai - a2)h\i^^y)p + 2(ai + a2)A\^,^y)p + -a4^''(^,^)p = 0, (17) 



and 



-p 4 4 

- (oi + a2)h [^,^]p (ai - -a^)nA^y - 2(03 + -a3)>l''[^,^]p 

-\aiA\,,^p + ^a^e^yxpd^d^i^^h'" - A^^) = 0, (18) 

respectively, in spatial regions far from the source. Here the d'Alembertian operator is 
given by □ = 9^*9^, hp,i, denotes 

= h^^u - ^Vi^uh, h = r]^,l,h'"', (19) 
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and Afj,^ stands for the dual of A^j^ defined by 

■^ixv — 2^ixv\p-^'^'^ ^ (20) 

where e^j^po- is the completely antisymmetric tensor normalized as e^^^"^ = +1. It is clear 
from (17) and (18) that the symmetric field h^i, and the skew-symmetric field A^,, are 
coupled to each other unless the parameters satisfy the condition ai + a2 = and 
a4 = 0. As is easily checked, the linearized field equations (17) and (18) are invariant 
under the gauge transformation 

h HI/ — hfj^n — 2^(^^,^), A ^j^i, — A^i, — (21) 

where the are small functions which leave the fields weak. By means of this freedom 
we can require the gauge condition 

= 0, (22) 
which we shall assume henceforth. Then Eqs. (17) and (18) become 

V - 2 (°i + a2){Vu,u^h - h^^u) + 2(ai + a2)A^(/.,,.)p + -a4A^(p,^)p = 0, (23) 

(ai - ^a^)uA^,, - 2(a2 + ^a3)^^[^,.]p - '^aiA\^^,^p - ^aief.^xpd^d^A"" = 0. (24) 

In order to calculate the total energy and momentum of the isolated system under 
consideration, it is enough to obtain asymptotic solutions of (23) and (24) up to order 
0{l/r). Multiplying on (23) [or (24)], we obtain 

□ [(ai + a2)A\ + a4Tj ^ = (25) 
with 04 = (04/3). Let us define B^j^^, — by 

Bf^v = (oi + 02)^,^1/ + claA^j,, (26) 
which can easily be solved with respect to A^j^^, to give 

^i^"" ^ (r. ^r. \2 ^7,2 [(^1 + «2)-B^i/ - a^^B^^ , (27) 
\ai -|- + ^4 '- -' 

under the assumption that (oi + 02)^ + a| 7^ 0, or equivalently (oi + 02) 7^ and/or 
a4 7^ 0. We note that when this assumption is violated, the theory is reduced to the 
special case of Hayashi and Nakano^^^ and M0ller.^'') 
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We now assume that the system under consideration is stationary and that the 
system as a whole is at rest somewhere around the origin. Then h^i, and Af^^, are 
time-independent. Rewriting (23) and (24) in terms of Bf^^,, we have after a shght 
modification 

Ah = 0, (28) 
^A^oo = 0, (29) 

^A/^Oa - 50/3,a/3 = 0, (30) 
/lA-Boa + f 2 Bop, (3a + h^al3-yBs[f3,-y]S = 0, (31) 

-Y^Ka + 2(^1 + «2)/i,a/3 + 25^(„,^)^ = 0, (32) 

fl^Ba/3 — '2f2Bs[a,l3]S + h^ap^B^s^S-y — 0, (33) 

where we have introduced /i, /2 and fz by 

4 

/i = (oi + a2)(ai - -03) + 04, (34) 
y 

4 

/2 = (ai + a2)(a2+ gOs), (35) 

4 

/a = 04(02 + gOa)- (36) 

We note that except for ASoq and AB^p Eqs. (30)~(33) involve Bqo, and S^/j only in 
the form of 3-dimensional divergence. 

According to (25) with (26) the divergences, -Boa,a and Bap,p, satisfy the Laplace 
equation, and therefore, they are given up to 0{l/r^) by 

Bap,p - ^2 ' y"^^) 

where B^ is a constant vector, and F^p is a nonsymmetric constant tensor. Here the 
radial unit vector is defined by n" = x"/r without making distinction between 
upper and lower indices. We have omitted from our consideration the possibility that 
the leading terms of the solutions (37) and (38) begin from 0(l/r). This is because if, 
on the contrary, they begin from 0(l/r)-terms, Eqs. (30)~(33) will have no solutions 
of h^i, and B^^, which tend to zero for large r. Since B^p is by definition antisymmetric 
with respect to a and /3, the constant tensor F^p must be of the form 

hp = F[ap] + ^F5^p (39) 

with F being the trace. With the help of (37)~(39) we can solve the field equations 
(28)~(33) by following the ordinary procedure to solve the Laplace equation. (See the 
Appendix for a summary of the procedure.) 
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^From (28) we see that h can be expressed by h = p/r up to 0(l/r), where p is a 
constant. Using this expression for h and Eq. (38) in (32), we obtain 

haP = + |4(55(a-F/3)7 - -FSa^S^S + P{ai + ^2) (SajS^s - -^^apS^S^ | 

' _ " (40) 

where Gai3 is a symmetric, constant tensor. Applying the condition (22) to ha/s of (40), 
we see that Gaj3 is given by 



2 



4 s s 1 



(41) 



Since Gap is symmetric, Eq. (41) imphes 

F^ap] = 0. (42) 

As for the field equation for we note that the second term on the left-hand 
side of (33) vanishes owing to (39) and (42). Applying the same procedure to (33), we 
obtain 

i?a/3 = ^ + ^^(^7 - 3n,n,5,), (43) 
r 6/1 r 

where we have assumed that the constant /i of (34) is nonvanishing. Here Hafi is an 
antisymmetric constant tensor to be fixed by the relation (38): We have 

Ha/3 = ^0/37-87, (44) 

F = 0, (45) 

the latter of which together with (39) and (42) implies that Bai3,/3 vanishes up to 
0(l/r2). 

We can analyze the field equations for hoa and Boa in the same manner. We here 
summarize the asymptotic solution up to 0(l/r) for hj^i, and Bj^i, as follows: 

h = l, (46) 

5» = ^^^1^, (47) 

hoa = + (48) 

— (ai + a2)p ,^ , , 

han = — — {da/3 + nan/3), (49) 
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Boa — 



B. 



a/3 



Ba + ^ (Sa + rianpBp) 



2/1 r 



(50) 



(51) 



where p is an unknown constant and Ba is an unknown constant vector. From (19) 
and (27), we obtain the most general expression for h^, and A^i, up to 0(l/r) 



h = — , 
r 



(52) 



^00 — 



3a2P 



hoa — 



2r 



2r ' 

p{ai + 02) {dai3 + nccUp) 



6ai 



2/ir 



2Ba ( «i - j + f a2 + ^aa j (b^ + naTipB^J 



A - . R 
fir 



(53) 
(54) 
(55) 
(56) 
(57) 



The covariant components of the parallel vector fields are then represented up to 0(l/r) 

by 

(58) 



b<\ = 1 - ^, 

r 



/i 



3ai 



/i 



2r 



(59) 



/i 



+ < 



3ai^ A I 2r 



m {(02 - 2ai)(5aa + (oi + a2)nana} , a^eaa^B^ 



9aia2r 



+ 



fir 



with m being a constant given by 



(60) 



(61) 



m 



3a2P 



(62) 



Thus the parallel vector fields are asymptotically characterized by the constant vector 

Ba besides the constant m at large spatial distances from an isolated, stationary system. 
The constant m can be interpreted as the gravitational mass of the source, as seen from 
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(53). The physical meaning of the constant Ba will become clear when we calculate in 
the next section the total energy and momentum of the isolated system. 

The parallel vector fields of (58)~(61) are not spherically symmetric because they 
involve the constant vector B^- This is to be contrasted with the metric in general 
relativity, which is governed by the Einstein equation. More precisely, if one selects 
a system of coordinates with respect to which the isolated system whole at 

rest, and if one neglects gravitational radiation by the system, the Einstein equation 
requires that the metric is spherically symmetric and static up to 0(1/ r) for large r. 
We note, however, that when Ba — 0, the parallel vector fields (58)~(61) reduce to 
the static, spherically symmetric solution expanded up to 0{l/r). 

We have been assuming that the parameters satisfy (ai + 02)^ + 04 7^ and /i 7^ 0. 
Now that we have obtained the solution of {b''^) up to 0(l/r), we can study the limiting 
cases a4 = 7^ (oi + 02) and (ai + 02) = 7^ 04. In both cases we can use A^i, directly 
instead of B^j^^. In the case 04 = and (ai + 02) 7^ 0, we can further take the limit 
(ai + 02) — > 0: It is easy to check that the resultant h^^ and A^i, satisfy Eqs. (23) and 
(24) with (01 + 02) = and 04 = 0. In case (01 + 02) = and 04 7^ 0, on the other hand, 
it is possible to take limit 04 — > only when (oi — 403/9) — > is taken simultaneously. 
The resultant theory is just the tetrad formulation of general relativity. 



|4. Calculation of energy and momentum 



Since the total Lagrangian is a scalar under general coordinate transformations, the 
total energy-momentum complex {TiZ+tf/") is represented as 

V^iT," + - W/^^, (63) 

with the superpotential U^'^^ being antisymmetric in u and A. Here t^^" is the canonical 
energy-momentum complex for the gravitational field derived from the gravitational 
Lagrangian density (4).^^^ The superpotential Un''^ is given by 




(64) 



which satisfies the M0ller condition (C). Since we have solved the field equation approx- 
imately at far distances, it is better to rewrite the superpotential (64) in the asymptotic 
form 
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K 



2 (ai - f ) ^m''^'" - («i - 2a2) (^^''^a^'^", - h''^\j^^,) 
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, (65) 



keeping only 0(l/r^)-terms. Here we use only Greek indices and raise (or lower) them 
with the Minkowski metric: Thus, for example, 6'^^ = S-^kb'^n and b//''^ — rj^'^dpbn''. he 
formula for the energy and spatial momentum are given by^^ 

= - lim / Uo^"dSa, (66) 

J r=const 

= lim / Ua'^dSf,. (67) 

J r=const 

The gravitational Lagrangian is assumed to be a quadratic invariant of the torsion 
tensor, and accordingly the canonical complex t^" behaves as 0(l/r'^) at far distances 
since the parallel vector fields are asymptotically given by (58)~(61). Thus, the M0ller 
condition (B) is also satisfied, and the energy and momentum are transformed as a 
4- vector. 

Now let us calculate the energy by using the asymptotic form of the parallel vector 
fields. The necessary components of the super potential are given by 



f . ^. (2^2 - aif\ 3aia4 ^ 

(oi + a2) H > mua + „ eapjUpB^ 

yaia2 I ^/i 



' Cn. - 1 , (68) 



where we have used the Newtonian approximation condition (10) to rewrite the first 
term inside the square bracket. It is easily shown that Eq. (68) satisfies Uq^"^^ = up 
to 0{l/r^) as it should. Using (68) in (66) we find that the last term of (68) does not 
contribute to the integral and that the energy is given by 

E^m. (69) 

This shows that the most general Lagrangian (4) including a parity-violating term is 
consistent with the equivalence principle. 

Next, we turn to the spatial momentum P^. The necessary components of the 
superpotential are given by 
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+2 (ai + (oi + a2)n[aB,3] - -^fimnpea/Sp - al (niaB/s] - Sa/sn^B^) 



(70) 



Taking the divergence of (70), we find that Ua^f^ ,p = 0, which means that the leading 
term of the Ud^^ p is of 0(l/r^), as it should. Now we are ready to calculate the spatial 
momentum to determine if it is vanishing or not. Using (70) in (67), we find after a 
lengthy calculation 

Pa = ^„ (71) 

which shows that the constant vector Ba has the physical meaning of the spatial mo- 
mentum of the isolated system. 

Since the isolated system is assumed to be stationary and as a whole at rest near 
the origin, its spatial momentum should be vanishing, and hence the constant vector 
Ba must be zero. This means that the divergence -Boa,a should be at most of 0(l/r^) 
for parallel vector fields at distances far from a stationary, isolated system. As we 
have seen in the previous section, however, such an asymptotic behavior of -Boq,^ does 
not follow from the linearized field equation in vacuum. Accordingly the asymptotic 
condition, 

limSo„,„ = 0(l/r^), (72) 

r— >oo ! V ' / V / 

must be imposed by hand, in addition to the condition (12) (or (13)). By contrast, as we 
noted below (45), the hnearized field equation implies the condition Bap^p — 0{l/r^). 

The asymptotic condition (72) is indeed satisfied when the gravitational field is 
weak everywhere and the linearized field equation, which is just given by Eqs.(17) 

and (18) with the right-hand sides being replaced by kT(^j,) and kTj^,^], respectively, is 
valid throughout whole space. Here T^^ is the energy-momentum tensor in the special 
relativistic limit and satisfies the conservation law di,T^^ = 0. Multiplying both sides 
of the linearized field equation by d'', and using the conservation law, we obtain the 
inhomogeneous Laplace equation for Boa,a 

— 2^0ap,aP, (73) 

where S^'^^ is the intrinsic spin tensor of the source, and we have used the Tetrode 
formula T^^'^ — {l/2)dxS^'^^ in the last step. For an isolated system confined to a finite 
spatial region, the solution Boa,a of (73) is of 0{l/r^) for large r. 

In the parity-conserving case with 04 = it is more convenient to use Aqq and A^p 
directly. Writing the asymptotic form of Aoa.aj which satisfies the Laplace equation, 
as BaUa/r'^, Eq. (71) reads 

Po, = (ai + a2)Ba. (74) 
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As explained at the end of §3, the solution of /i^i, and A^^, has a well-defined limit 
when we put — and then (ai + 02) = 0. Although the resultant solution involves 
the constant vector Ba, the total spatial momentum is vanishing as is shown in (74). 
This situation can be understood as follows. In the special case (ai + 02) = = 04, 
the linearized field equation for A^,^ is decoupled from that for h^i,, and therefore, the 
former equation is invariant under the gauge transformation 

^ fj,u = Afj,!/ + C[^i,l'] (75) 

with representing arbitrary small functions. We can use this freedom to choose a 
gauge in which Aoa and Aap are vanishing up to 0(l/r). Thus, in this special case, the 
constant vector Bg, is unphysical and can be eliminated by a gauge transformation. 



§5. Conclusion and discussion 



We have studied the energy and momentum of an isolated system in the tetrad 
theory of gravitation, starting from the most general Lagrangian that is quadratic in 
torsion and involves four unknown parameters ai, 02, 03 and 04, the last of which is 
associated with a parity-violating term. As a first application we considered the static, 
spherically symmetric case, where the parallel vector fields take a diagonal form. The 
solution of the field equation in vacuum is found to be the same as the exact solution of 
the parity-conserving case with 04 = 0.^^^ This is due to the fact that the axial- vector 
part of the torsion tensor is identically vanishing for diagonal parallel vector fields. 

The total energy and momentum of an isolated system are expressed by a surface 
integral over a large closed surface enclosing the system. It is then sufficient to know 
the asymptotic form of the parallel vector fields at distances far from the source, where 
the gravitational field is weak. In view of this, we analyzed the linearized field equation 
in vacuum which follows from the most general gravitational Lagrangian. 

It is well known that in general relativity the Einstein equation ensures that the 
metric tensor at far distances is the same as the Schwarzschild metric up to 0(l/r) 
for any isolated stationary system. This is usually shown^^^ by solving the linearized 
Einstein equation in vacuum, which is the Laplace equation in stationary case and can 
easily be solved. 

In our case, the linearized field equation in vacuum consists of 16 equations for the 
symmetric field /i^i^ and the skew-symmetric field A^i,, and is invariant under a gauge 
transformation, which allows us to impose the harmonic condition on the h^^. It is 
found, however, that the symmetric part and skew-symmetric part of the linearized 
field equation are coupled with each other unless the parameters satisfy Oi -|- 02 = = 
04. Nevertheless, we can solve the coupled equation up to 0(1/ r) assuming that the 
parameters satisfy (ai -|- 02)^ + (04/3)^ 7^ and (ai + a2)(ai — 4a3/9) -|- (04/3)^ 7^ 0. It 
is found that the general solution has two constants, the gravitational mass m and a 
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constant vector Ba- The general solution of the linearized field equation up to 0(l/r) 
is, therefore, not spherically symmetric, due to the precense of Ba- 

Using the definition of energy-momentum complex given by M0ller,^^'* we derive 
the superpotential U//^ from the general Lagrangian. It transforms as a tensor density 
under general coordinate transformations. Using the general solution obtained, we 
give an explicit expression for the components Wo°" and Uj^'^, which are shown to be 
divergenceless up to 0(l/r^), implying that the leading term of the divergence of those 
components is of order 0{l/r'^). 

From the components Wo°" we calculate the total energy, and show that the result 
is equal to the gravitational mass of the isolated system, i.e. E — m. Using the 
components Uj^^ , we then calculate the spatial momentum, and find that Pa — Bq.- 
Thus, we arrive at the result that the general solution at far distances is characterized by 
the total energy and the spatial momentum of the isolated system under consideration. 

The asymptotic solution with vanishing Ba is acceptable as a description of the 
parallel vector fields far from a stationary isolated system which is at rest as a whole. 
It indicates that the total momentum of the isolated system at rest is vanishing. The 
asymptotic solution coincides with the exact solution with spherical symmetry up to 
0(l/r). 

It is not yet clear whether or not the asymptotic solution with Ba 0, actually 
describes the asymptotic behavior of exact solutions. Such exact solutions, if they exist, 
should describe exotic systems which have nonvanishing total momentum although at 
rest as a whole. 

The vanishing of the constant vector Ba implies that the divergence -Boa.a of -Boa = 
(oi + a2)AQa + (04/3)^00 behaves at most as 0(l/r^) for large r. This asymptotic 
behavior does not follow from the linearized field equation alone. By contrast, the 
asymptotic behavior of Ba/3 — (oi + 02)^0/3 + (04/3)^^/3 is more severely governed by 
the linearized field equation: In fact, it is shown that Ba/3,13 vanishes up to 0{l/r^). 

For a weakly gravitating source, for which the field is weak everywhere and the 
weak field approximation can be applied, we can show that B^a,a behaves as 0(l/r^). 
For the fully relativistic case we do not know the reason why Boa,a should behave in 
this manner. Therefore, further study is required to establish the asymptotic condition 
of the parallel vector fields, in particular, of the skew-symmetric field A^^. 

The study of the weak field case also needs more investigation. Although in the 
present paper we solved the linearized field equation for the stationary case, the same 
procedure may be applicable also for the time-dependent case. This will be studied in 
future work. 

The present analysis is based on the assumption that the linearized theory can be 
applied at distances far from the source. The solution obtained by Mikhail ct al.^^^ 
which violates E = m, does not satisfy this assumption. It can also be shown that 
their solution does not satisfy the M0ller condition (B), because the divergence Uj^'^ 
is 0(l/r^/^) for their solution. 
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Appendix A 

Solution of an Inhomogeneous Laplace Equation 

Here we derive the general solution of an inhomogeneous Laplace equation of the 
form 

1 oo 

AU{x,y,z) = - E A„,a....a;v^"'^"' • • (A • 1) 

^ N=l 

where Aaia2--ajv totally symmetric and traceless, and represents for the unit vector 
x°'/r. We look for a solution in the form 

-j^ oo 

U{x, y,z)^-J2 5aia....a^n"in'*^ • • • n"- , (A • 2) 

^ Ar=o 

where Baia2-aN totally symmetric and traceless. Applying the Laplacian operator 
on (A-2), and substituting the result in (A-1) we obtain 

U{x,y,z) = E T77T7— TT^aia^-a.n^^n"^ • • - n"-, (A • 3) 

r r N{N + 1) 
where B is an arbitrary constant. Here we have used the following proposition: 

when Aaiaj - ajv is totally symmetric and traceless. 
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